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Abstract. - We examine the impact of a solid sphere into a fine-grained granular bed. Using 
high-speed X-ray radiography we track both the motion of the sphere and local changes in the 
bed packing fraction. Varying the initial packing density as well as the ambient gas pressure, 
we find a complete reversal in the effect of interstitial gas on the impact response of the bed: 
The dynamic coupling between gas and grains allows for easier penetration in initially loose beds 
but impedes penetration in more densely packed beds. High-speed imaging of the local packing 
density shows that these seemingly incongruous effects have a common origin in the resistance to 
bed packing changes caused by interstitial air. 



Introduction. — Granular materials often exhibit be- 
havior intermediate between that of conventional solids 
and liquids. Probing the resulting combination of liquid- 
and solid-like properties a number of recent studies in- 
vestigated the impact of a large object into a bed of dry 
grains. These studies focused on issues such as the drag 
on the impacting object [T]^6] , crater formation [7[|lQ] the 
corona-like splash formed immediately after the impactor 
hits the bed surface [To], and the subsequent jet of grains 
formed by the collapse of the cavity left by the impactor 
So far however, 
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almost all work considered the 
limit of loosely packed, marginally stable beds that read- 
ily compact in response to perturbations. On the other 
hand, densely packed beds must dilate in order for grains 
to move out of the way of inserted objects. This implies 
different resistance not only for slow, quasi-static pertur- 
bations [17] but also suggests that there should be a sig- 
nificant change in the dynamics for faster impacts. 

An important feature of the impact dynamics in gran- 
ular systems is the coupling between the interstitial gas, 
typically air, and the grain packing. For fine grained beds 
(grain diameters below ~ 150 fim) this interaction can 
drastically change the impact dynamics. In particular, in 
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the presence of interstitial gas an impacting sphere pene- 
trates much deeper than in vacuum [l4l [T8]. O ne possible 
explanation, suggested by Caballero et al. |14|, is drag re- 
duction due to local fluidization provided by a layer of gas 
flowing around the moving sphere. A different explana- 
tion relies on the presence of interstitial gas throughout 
the bed. In situations where the gas permeability is low, 
i.e., fine-grained beds, the gas is effectively trapped during 
the short impact duration and can resist global packing 
density changes 
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For loose packings this leads to a 
response whereby a large portion of the bed surrounding 
the impact behaves similar to an incompressible liquid, al- 
lowing the impactor to sink in [18]. On the other hand, 
for denser packings the same mechanism predicts that the 
presence of gas should impede the impactor movement, 
since the gas-particle coupling should tend to counteract 
dilation. 

Here we test for these two different scenarios by varying 
both the interstitial air pressure and the initial packing 
of the bed. Using x-ray radiography, we track both the 
motion of the sphere through the bed and resulting local 
packing densities changes. 

Setup. — X-ray imaging was performed at the 
GeoSoilEnviroCARS beam line at the Advanced Photon 
Source as in refs. [l3][T5][T8] . The experimental setup and 
image processing are detailed in 15 , and we only out- 



p-1 



J. R. Royer et al 



line key aspects here. X-ray transmission through the bed 
was recorded at 6000 frames per second with resolution of 
29 /xm/pixel. The beam size restricted the field of view 
to 22 mm x 8.7 mm sections of the bed, so in order to 
capture the dynamics across the full vertical extent of the 
bed, movies of multiple, independent sphere drops, were 
aligned and synchronized. The detector was calibrated 
pixel by pixel to convert intensity to packing density, al- 
lowing us to correct spatial variations in beam intensity 
and detector sensitivity. 

For each experiment, a steel sphere (diameter D s = 
12mm) was dropped from a height of 35 cm into a bed 
of boron carbide (B4C) particles (diameter d = 50 /im 
± 10/im). The bed was contained in a cylindrical poly- 
carbonate tube with 35 mm inner diameter. Before each 
drop the bed was aerated by dry nitrogen entering through 
a diffuser built into the bottom of the container. After 
slowly turning off the nitrogen flow, the bed reproducibly 
settled to a loosely packed state with an initial average 
packing fraction O = 0.51 ±0.01. In this loose state the 
bed depth was 9 cm. For experiments with denser beds, we 
compacted the bed by gently tapping the chamber walls 
until the top surface of the bed fell to fixed level. The 
densest bed we obtained was 7.5 cm deep, corresponding 
to a 17% decrease in volume. With the x-rays we mea- 
sured the initial packing of the dense bed to be 0o = 0.60 
±0.01, in agreement with the change in height. In both 
the loose and compacted beds, the initial packing varied 
by no more than 1-2% across the bed height. The sys- 
tem could be sealed and evacuated down to pressures of 
0.7 kPa. The pump speed was limited to prevent air from 
bubbling up and disturbing the packing. 

Results. — In Fig. [l]we show the motion of the sphere 
and the local change in bed packing density A(j) = (j)(t) — (j)Q 
during penetration of initially loose and initially dense 
beds at atmospheric pressure (Pq = 101 kPa) and in vac- 
uum (Pq < 1 kPa). These images reveal striking differ- 
ences in the impact dynamics. 

In the initially loose bed under atmospheric pressure 
(Figs, [l] (a) - (e)) the sphere easily sinks into the bed, 
leaving a large cavity behind it. In fact, the sphere reaches 
the bottom of the 9 cm deep bed and even bounces before 
coming to rest. In the loose bed under vacuum (Fig. [lf-j) 
the sphere still penetrates into the bed and opens up a cav- 
ity, but the resistance of the bed is substantially increased. 
In contrast to the loose bed under atmospheric pressure, 
the sphere only sinks 3.5 cm below the top surface and the 
resultant cavity is much smaller. 

In the initially dense bed the resistance to penetration 
is significantly increased, stopping the sphere before it can 
even sink below the top surface. However, now the effect 
of the interstitial gas is reversed: in vacuum the sphere 
sinks nearly twice as deep (Fig. [l]>t) as in atmospheric 
pressure (Fig. [ljt-o). To quantify this pressure depen- 
dence, we track the sphere position z s (t) varying the am- 
bient pressure (Fig. [3] a). These trajectories show that the 
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Fig. 1: (Color) Composite x-ray images showing changes 
in local packing fraction Acj) = <j) — 0o after a metal sphere 
impacts a granular bed surface at time t — 0. The im- 
pact sequences (left to right) contrast initially loose beds, 
(a)-(e) and (f)-(j), with initially dense beds, (k) - (o) and 
(p) - (t). For each the response under atmospheric pres- 
sure (lOlkPa) as well as vacuum (lkPa) are shown. The 
impacting sphere appears black while the cavity behind 
the sphere appears white. Supplemental movies for each 
impact sequence are available online: Loose_atmospheric.mov 
(Quicktime 2.1 MB), Loose_vacuum.mov (Quicktime 2.3 
MB), Dense_atmospheric.mov (Quicktime 852KB) and 
Dense_vacuum.mov (Quicktime 373KB). 
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Fig. 2: Transition in impact dynamics. The final penetration 
depth of the sphere, Zf, as a function of initial packing fraction 
0o at Po = 101 kPa (solid symbols) and in vacuum Po = 0.9 
kPa (open symbols). For the loose (0o = 0.51) bed at 101 kPa, 
the sphere hit the bottom of the bed in the chamber used for x- 
ray measurements so z s f was measured separately in a deeper 
bed using a line attached to the end of the sphere. Inset: The 
same data on a log scale to highlight the crossing around 0o = 
0.58. 
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Fig. 3: Sphere dynamics. (Color Online), (a) Vertical position 
z s (£) of the bottom tip of the sphere as a function of time after 
impact (t = s) for an initially dense bed 0o = 0.61. From 
top to bottom: P = 101 kPa, 50.6 kPa, 14.3 kPa, 6.7 kPa, 3.6 
kPa and 0.7 kPa. Inset: z s (t) for an initially loose bed 0o = 
0.51 and from bottom to top: P = 101 kPa, 12 kPa, 8.7 kPa, 
4.9 kPa and 0.7 kPa. (b) Velocity v s (t) computed from curves 
in (a) for 0o = 0.61. Inset: Comparison of v s (t) in vacuum 
(Po < 1 kPa) for initially loose and compacted beds. 
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Fig. 4: Packing changes in vacuum, (a-d) Space-time plots 
of the centerline of the x-ray movies in Fig. 1, showing the 
change in packing A0 in vacuum (Po < 1 kPa) for different 
initial packing densities (a) 0o = 0.51 =b 0.01(b) 0o = 0.53 =b 
0.01(c) 0o = 0.58 ± 0.02, (d) 0o = 0.60 ± 0.01. (e-f) Packing 
density observed below the sphere after it has come to rest: 
peak change (e) and peak value (f) for different initial packing 
fractions 

penetration depth monotonically increases with decreas- 
ing ambient pressure Po in dense beds. In loose beds, by 
contrast, the penetration depth increases with increasing 
ambient pressure (Fig. [3] a inset) [l4{[l8] . 

To further examine this reversal in the impact dynamics, 
we plot the final sphere depth Zf against 0o for beds at at- 
mospheric pressure and under vacuum in Fig. |2j At atmo- 
spheric pressure, the final sphere depth depends strongly 
on the initial packing density, decreasing by over an or- 
der of magnitude as this density is increased from 0o = 
0.51 to 0.62. Under vacuum, on the other hand, the pen- 
etration depth is nearly independent of </>o, with only a 
slight increase at our lowest value 0o = 0.51. The inset in 
Fig. [2] shows the same data on a log- log scale to highlight 
the crossover from deeper penetration under atmospheric 
pressure to deeper penetration under vacuum. Though 
our bed preparation method is not optimal for controlling 
the packing fraction precisely and we only have results 
for four values of 0o, the data indicate that the crossover 
occurs near O = 0.57-0.58. 

In the initially loose bed, during impact under vacuum a 
large front of compacted grains builds up ahead the sphere 
(Fig. [lf-j). The extent of the compacted region grows as 
the sphere moves through the bed, so the edge of this 
compaction front travels faster than the sphere. However, 
after the sphere comes to rest the front does not continue 
to propagate and instead stops before reaching the bottom 
of the bed (Fig. [l]j). This differentiates this compaction 
front from propagating pulses observed, e.g., in granular 
chains [l9]. The maximum packing density in the com- 
paction front is ~ 0.54, much less than the packing 
densities of up to ~ 0.61 obtained by manually tapping 
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the chamber to compact the bed. 

In the initially dense bed under vacuum, there is con- 
siderable dilation around the moving sphere, as one would 
expect [2Q|[21] . Just like the compaction front in the loose 
bed, the dilation front does not propagate through the bed 
after the sphere has come to rest. The packing fraction 
just below the sphere in the dilated region is about <p ~ 
0.56, just slightly higher than the 0.54 measured for the 
compaction front in the loose bed. 

For both loose and dense beds under atmospheric pres- 
sure, most of the changes in packing density occur in a 
small region just ahead of the sphere (Figs. [T^i-e, k-o). 
This is in line with the idea that the interstitial air re- 
sists changes in the packing density. An initially loose, 
marginally stable bed thus stays liquid-like and flows out 
of the way of the sphere [18], while an initially dense bed 
stays solid-like and flows very little. 

To quantify how the resistance of the dense bed to pen- 
etration depends on the ambient pressure, we compute 
v s = dz s /dt from the measured trajectories z s (t), shown 
m Fig. |3]b. There is a sharp, initial deceleration which is 
nearly independent of the ambient pressure. At reduced 
ambient pressure, at later times there is turnover to a 
lower deceleration. In the inset in Fig. ||b we compare 
the sphere velocity in vacuum in an initially loose and ini- 
tially compacted bed. After the initial rapid deceleration 
in the dense bed, the deceleration at later times is very 
close to the deceleration in the loose bed. 

The connection between packing density changes and 
impactor motion is most clearly seen in space time plots 
from composite x-ray movies under vacuum at varying ini- 
tial packing fractions (Fig. |4^i-d). In the loosest bed Fig. 
|4^i there is a large, clear compaction front which extends 
far ahead of the sphere, but this front dies out and stops 
at roughly the same time the sphere comes to rest. As the 
initial packing is increased (Fig. [4j)), the sphere does not 
penetrate as deep and the compaction front becomes less 
pronounced until around 0o =0.58 there is almost no local 
density change as the sphere move through the bed (Fig. 
[4]:). Finally, in the densest bed ( Fig. [4]i) we see that the 
dilation front ahead of the sphere stops the same time the 
sphere comes to rest. Plotting the peak change in packing 
below the sphere (Fig. |4^) we can track this change from 
compaction in loose beds to dilation in dense beds. At the 
same time, the magnitude of the packing fraction at this 
peak does not exhibit any clear dependence on (/>o, varying 
between 0.54 and 0.58 (Fig. g). 

Discussion and Conclusions. — Our results demon- 
strate that, during sudden impact, the primary role of 
the interstitial air in a fine-grained granular bed is to op- 
pose changes in the bed packing density. In loose beds 
the trapped gas prevents compaction of the grain arrange- 
ment, and as a result the bed flows out of the way of the 
intruder, similar to an incompressible fluid. This allows 
for relatively easy penetration and explain why a heavy 
sphere, even if started just above the bed with zero impact 



velocity, can sink in deeply 12 . In sufficiently dense beds, 
the trapped gas has the reverse effect, working against di- 
lation of the grain arrangement and thereby resisting pen- 
etration. 

In the absence of interstitial gas, the impacting object 
can change the grain packing density much more effec- 
tively, generating a pronounced compaction front if the 
bed is loose and dilating the bed when it is dense. For 
loose beds, in particular, the result is behavior similar to 
that observed ahead of a plow pushed into soil or snow, as 
seen by the reddish-colored regions of enhanced density in 
Figs, lf-j and 4a. 

If local, gas-mediated fluidization around the impact- 
ing sphere were to provide the mechanism for deep pen- 
etration in loosely packed beds, then the presence of gas 
should have allowed for deeper impact also in more densely 
packed beds. The fact that we find the opposite behav- 
ior, whereby gas reduces the impact depth in the densely 
packed case, appears to rule out this scenario. Indeed, in 
qualitative terms the same reversal in relative penetration 
depth with increasing packing density was already pointed 
out in a NASA Technical Note in 1963 [22] 

As the packing density is increased, we observe a grad- 
ual transition from compacting to dilating response of the 
bed. For the type of bed material examined here, the 
crossover between these two regimes occurs at a packing 
fraction <\> ~ 0.58 (Fig. 4e). Within our experimental 
accuracy, this agrees well with the packing fraction <j> ~ 
0.57-0.58 at which the penetration depth with ambient 
air present becomes smaller than the same depth mea- 
sured in vacuum (Fig. 2 inset). This finding says that the 
crossover from compaction to dilation occurs at, or at least 
very close to, the density at which the role of interstitial 
gas is negligible. It suggests, therefore, that the partic- 
ular crossover value should be independent of gas-grain 
interactions. 

In principle, this makes comparison possible with exper- 
iments or calculations where the coupling to air does not 
need to be considered, such as slow, quasi-static penetra- 
tion. On general grounds we expect the crossover density 
to lie above the value (j) r i p for random loose packing, since 
(j) r i p corresponds to marginally stable configurations that 
collapse and compact when perturbed 23 . For slowly 
settled non-cohesive spheres, <j) rcp ~0. 56-0. 5 7, with lower 
values for rougher particles and larger values for larger 
mismatch with the suspending fluid 21p3 . A few percent 



above (j) rcp Schroter et al. indeed observed a transition in 
the resistance of a granular bed to quasi-static penetra- 
tion by a rod 17 . Using spherical glass beads Schroter 
at al. found this transition to occur at values just below 
0.60. In our experiments, the deviation from sphericity 
in the boron carbide particles is likely to be responsible 
for pushing <j) r i p to lower values and with it the crossover 
density to around 0.58. 
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